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Abstract 

We prove a regularity result for the unstable elliptic free boundary problem 

Am = -X{«>0} (0.1) 

related to traveling waves in a problem arising in solid combustion. 
The maximal solution and every local minimizer of the energy are regular, that is, {u = 0} 
is locally an analytic surface and locally analytic functions. More- 

over we prove a partial regularity result for solutions that are non-degenerate of second 
order: here {u — 0} is analytic up to a closed set of Hausdorff dimension n — 2. We 
discuss possible singularities. 
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1 Introduction 



We prove a regularity result for the unstable elliptic free boundary problem 

Au = -X{„>o} ■ (1-2) 

The problem (|1.2|) is related to traveling wave solutions in solid combustion with ignition 
temperature: let us consider the solid combustion system 

d t 9 - A9 ^ g(i 1 )x{e>e }, , ± ^ 

dtv = g(v)X{e>e } 

where 9 is the temperature, 9q is the ignition temperature and r\ S (0, 1) is the fractional 
conversion. Although an ignition temperature has no meaning for gas flames, it has 
been recently rediscovered and used in combustion synthesis (see for example |20].|3|.p] 
and HZP- The reaction kinetics suggested in 1201 p. 1462] is 51(77) = k X{ t? <i}i but the 
particular form of g is not of importance for what follows. Solving the ODE we obtain 
for u := 6 — 9q 

d t u -Au = cx{ u>0 } , (1.4) 

where c is a non-negative memory term depending on (t, x) as well as the history 
{u(s,x),s < t}. Although c is important when considering the large-time behavior of 
u, we may consider <|1.4ll to be a perturbation of 

d t u -Au = X{u>o} (1-5) 

when we are interested in transient or local phenomena at the ignition front d{u < 0}. 
Actually the traveling pulses in our model correspond well to the fingering phenomenon 
for burned regions observed in solid combustion experiments (see for example |24| ) . 
We are interested in traveling wave solutions. As we deal in the present paper mostly 
with regularity issues, we may drop the drift term resulting from the time derivative of 
the traveling wave. 

An equation similar to our elliptic one arises in the composite membrane problem (see 
jl(J|. jSj). Another application is the shape of self-gravitating rotating fluids describ- 
ing stars (see equation (1.26)]). 

From a mathematical point of view, (|1.5|l is the equation of the parabolic obstacle prob- 
lem with inverted sign, and (|1.2f> is the equation of the elliptic obstacle problem with 
inverted sign. The change of sign changes the character of the problem drastically in that 
it changes the stable obstacle problem into an unstable problem. In l|1.2f> and l|1.5l) . we 
find examples of non- uniqueness, bifurcation phenomena etc. 

As surprisingly many known free boundary problems turn out to be stable problems, 
this means, unfortunately, that many known methods in free boundary problems do not 
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apply here. Examples of PDE techniques that do not work are, apart from all one-phase 
methods, the Bernstein technique, the Alt-Caffarelli-Friedman monotonicity formula ([3]) 
and the differential inequality technique of Cazenave-Lions ( 8 ). 

Our main result is that the maximal solution and every local minimizer of the energy 
are regular, that is, {u = 0} is locally an analytic surface and u\ ^ u>Q y , u \ ^ u<Q y are locally 
analytic functions (Theorem |OJ. 

The surprise is that - in contrast to the usual procedure - we obtain C 1,1 -regularity of 
local minimizers by proving regularity of the free boundary first! 

For general solutions that are non-degenerate of second order, we prove a partial regular- 
ity result: here {u = 0} is smooth up to a set of Hausdorff dimension n — 2 (Proposition 
16. 311 . We discuss the behavior at possible singularities. 

In case of a non-degenerate minimal solution in two dimensions we also obtain that 
{u = 0} consists of Lipschitz arcs meeting at right angles in at most finitely many sin- 
gularities. 

Acknowledgment: We thank Ivan Blank, Carlos Kenig, Herbert Koch and Stephan 
Luckhaus for fruitful discussions. 

2 Notation 

Throughout this article R" will be equipped with the Euclidean inner product x ■ y and 
the induced norm \x\ . We define ei as the i-ih unit vector in R™ , and B r (x°) will denote 
the open n-dimensional ball of center a; , radius r and volume r n cu n . We shall often use 
abbreviations for inverse images like {u > 0} :— {x G il : u(x) > 0} , {x n > 0} := {x G 
R n : x n > 0} etc. and occasionally we shall employ the decomposition x = (xi, . . . , x n ) 
of a vector x G R™ . We will use the fc-dimensional Hausdorff measure 7i k approximated 
by Ti, k ' S which we define as the of |14j . When considering a set A , \A shall stand for 
the characteristic function of A , while v shall typically denote the outward normal to a 
given boundary. 

3 Existence and Non-Degeneracy 

Lemma 3.1 (Existence of a maximal and a minimal solution) 

Let f2 be a bounded domain of class C 2 ' a in R", and assume that the Dirichlet boundary 
data u D G C 2 ' a (n). 

Then there exists a maximal solution u with the following properties: u £ W 2,p (fl) for 
every p G [1, +oo), u = ur> on dfl and u > v for every subsolution v G W 2 ' n (fl) of 
in Q' C f2 satisfying v < u on <9f2'. 

There also exists a minimal solution with analogous properties. 
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Proof: We prove the existence of the maximal solution: We consider a regularization of 
the equation from above, 

Av = -AO), 

where /3 e £ C°°(R), (3 e (z) > X{z>o} —'■ P( z ) m R and /3 e I P as e j 0. By Perron's method, 
there exists a maximal solution u e in the above sense for each e > 0. By W 2 ' p -estimates 
the family (ii e ) ee ( ,i) is bounded in W 2 ' P (Q). Moreover, for any subsolution v of 11.2fl 
and < £ < £ we obtain 

Av > -f3(v) > -j3 c {v) and Au~ e = -fe{u t ) > -/? e (ttg) 

so that v and ug are subsolutions of the e-equation. Consequently v < u e and ui < u e . 
Using the fact that D 2 u — a.e. on {u = 0}, it follows that u(x) := lim e ^o u t (x) is the 
maximal solution. 

Lemma 3.2 There exists a positive constant c n depending only on the dimension n, such 
that the maximal solution u with respect to given boundary data satisfies 

B r (x°) C tt, inf u > ~c„r 2 =>• u(x°) > . 

dB r (x°) 

Thus for any x° £ {u = 0} and r < dist(ir , Q c ) we obtain ii)igB r (x ) u < ~c n r 2 . 
Proof: Compare u(x° + rx)/r 2 for suitable to the "stationary pulse" p(9x)/0 2 , where 



p(x) 

in the case n — 2 and to 

p(x) = ■ 

in the case n > 2. 



(l-\x\ 2 )/4, \x\<l, 
-log(M 2 )/4, \x\>l 



^(1-N 2 ), \x\<l, 
^(\x\ 2 --l), \x\>l 



Remark 3.3 For the case of two dimensions the optimal constant ci = l/(4e). 
Definition 3.4 (Non-degeneracy) Let u be a solution of in fl, satisfying at x° £ 

n 

liminf r- 2 ( [ u 2 dH^ 1 ) > . (3.6) 

\JdB rm ( X °) J 

Then we call u "non-degenerate of second order at x . We call u "non-degenerate of 
second order" if it is non- degenerate of second order at each point in tt. 
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Remark 3.5 The maximal solution is non- degenerate of second order. 



Lemma 3.6 Each minimizer u of the energy 

E(v) = / |Vw| 2 - 2max(v,0) 

JB ro (xO) 

in K := {v £ W ' (B ro {x )) : v — u on dB rQ (x )} satisfies at x° the second-order 
non- degeneracy property 

Proof: Define for r £ (0,ro) the solution v := u(x° + rx)/r 2 and let p be the stationary 
pulse in £?i(0) with boundary data infg Bl ( ) v. Comparing the energy of v to that of 
in := max(p, v) we obtain 

0> / |Vu| 2 - |Vw| 2 - 2max(i>,0) + 2max(w,0) 

•/-Bi(O) 

-(Au + Aw)(v — w) + 2(max(w, 0) — max(«, 0)) 

JJyJU 

V + p . 



's 1 (o)n{w<o}n{p>o} 

Now assume that v m is a sequence of minimizers such that v m (0) = and mfgs^o) v m — > 
as m — > oo. Let e e (0, 1) be fixed. In the case that £ n (.Bi_ £ (0) n {v m < 0}) ■/* 
as m — > oo, we obtain immediately a contradiction since p m is for large m on _Bi_ e (0) 
uniformly estimated from below by a positive constant depending only on n and e. In the 
case £™(£?i_ e (0) n {v m < 0}) — > as m — > oo, we know by lin^^oo inf Sl ( ) v m — > and 
IdB (o) Wm dTt" 1 " 1 < for every s € (0, 1] (which follows from v m being superharmonic 
and at the origin) that J* Sl ( ) -> as oo. Consequently 



0<— / Au m = / ~X{v m >o} < -ci < as m — ► oo, 

JBi_,(0) JSi_ e (0) 

a contradiction. Thus w(a; ) = implies 

inf -i — ^ i < -c 2 < 

Bi(o) r z 

for all r £ (0, ro). 

Lemma 3.7 The maximal solution u ma x is the maximal minimizer of the energy 

E(v) = / |Vu| 2 - 2max(v,0) 
Jn 

in K :={v £ W h2 (Q) :v = u D on dtt}. 
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Proof: For any v G K we get 



> / max(ti,0) — max(« 
Jn 




0) + WmaxX{„>0} ~ V X{ 



"max>0} 



(X{i»o} + X{« max >o})("max -v) + 2(max(i;, 0) - max(u max , 0)) 



n 




n 



|Vu max | 2 - |Vu| 2 +2max(t),0) - 2max(u max ,0) . 



4 Monotonicity Formula and Frequency Lemma 

A powerful tool is now a monotonicity formula introduced in |22| by one of the authors 
for a class of semilinear free boundary problems. For the sake of completeness let us 
state the unstable case here: 

Theorem 4.1 (Monotonicity formula) Suppose that u is a solution of fl.ty) in f2 and 

that Bs(x°) C fl . Then for all < p < a < S the function 



We also need the following frequency lemma, which has been proven in |211 Lemma 4.1]: 

Lemma 4.2 (Frequency lemma) Let a-lgN, let w £ W 1 ' 2 (Bi(0)) be a harmonic 
function in -Bi(O) and assume that D J w(0) = for < j < a — 1 . 



and equality implies that w is homogeneous of degree a in B\ (0) . 

5 Classification of Blow-up Limits 

A result related to the following classification of blow-up limits is contained in jSJ Theorem 
2.5]. Note however that the proofs of the related parts are largely different. 





defined in (0, 5) , satisfies the monotonicity formula 
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Proposition 5.1 (Classification of blow-up limits with fixed center) Let u be a 

solution of in il and let us consider a point x a € f2 H {u = 0} R {Vit = 0}. 

1) In the case $ x o(0+) = — oo, lim r ^ r _3 ~" f gB , -.u 2 dH n ~ 1 = +oo, and for S (x° , r) := 

u(x° + rx) 
S{x°,r) 

as r — > is a homogeneous harmonic polynomial of degree 2. 

2) In the case $3.0(0+) £ (-00, 0), 

u(x° + rx) 

U ^ X ) : = ~2 

is bounded in W ' (Bi(0)), and each limit as r — > is a homogeneous solution of degree 
2. 

3) Else $ x o(0+) = 0, and 
u(x° +rx) 12 



in W^iBtiO)) asr^O 



Proof: In all three cases 
1 

'2 



d r \ [ ufdH"' 1 = [ UrdrUrdH"- 1 

= — I u r (\7u r ■ x — 2u r ) dTC 1 ^ 1 
r Jasi(o) 

- ( / -max(u r ,0) + |Vu r | 2 - 2 f u 2 r dU r ' 

r \JB 1 {0) JdB^O) 

1 

r 



$x°( r ) + / max(u r ,0) 

'Bi(O) , 



In particular, 



/ max(« r , 0) > M := -$ x o(0+) implies d r / u 2 r dW 

JBi(O) JdBxiO) 



> . (5.7) 



Observe now that 



f max(-u r , 0) 2 dK 1 - 1 < G x f 1 + / max(n r , 0) 2 dW n_ 



(5.8) 



supposing towards a contradiction that this is not true, f dBl r ) w 2 d7i n ~ l must be un- 
bounded for a sequence (f*ra) me ]\j not satisfying 1)5. 8fl . We may divide the function of 
the monotonicity formula 14. II 

/ |Vur m | 2 -2 / u 2 m cM"- 1 < $ a o(l) + 2 / max(u rm ,0) 
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by / aSl ( ) max(— u Vm , 0) 2 dH n 1 . As a subsequence r m — > 0, we obtain a weak L 2 (<9Bi(0))- 
limit u of i> m := u Tm j I Jq Bi ^ max(— u Tm , 0) 2 dTi. 71 ^ 1 J such that by |^r m | < 

^ 2 (jaBi(o) \ Ur rJ\ 2 dH n ^ 1 ^j (which follows from min(u rm ,0) 2 being subharmonic and 
u Tm being superharmonic and at the origin) 

|Vt; m | 2 -2 / v 2 m dH n ~ x < o(l) as m -4 00 , (5.9) 

Bl(0) JdBi(0) 

and v is in -Bi(O) a non-positive harmonic function satisfying v(0) — 0, implying by the 
strong maximum principle that v = in Bi(0). That poses a contradiction to the (by 
strong L 2 (<9-Bi(0))-convergence of w m and the fact that 

liminf / v^dH 71 - 1 > . 
Note also that in the case $> x o(0+) > —00, for small r and f g (r/2,r) 

1 > $ x o(r) - $ x o(f) = ;:2 S / 9Bi(0) (9 sUs ) 2 dW"- 1 

— /9Bi(0)( Ur ~ Uf) 2 dH n ~ 1 . 
Combining 15.10fl . (|5.8() and l|5.7(l we see that in the case $3,0(0+) > —00, 

/ u 2 > M implies d r [ u 2 dU^ 1 > . (5.11) 

But then Jg Bl m u 2 - dTC 1 ' 1 has to be bounded in the case $3.0(0+) > -co. 

From the monotonicity formula Theorem 14. II we infer that in the case $3.0(0+) > —00, 

each limit uq of u r in -Bi(O) is a homogeneous solution of degree 2. 

Moreover, in the case $3,0(0+) > 0, we obtain that each limit uq satisfies 

0<$ x o(0+) = f |V Uo | 2 -2max( Uo ,0) - 2 f u 2 dU n - 1 = - / X{u >o} ■ 

JBi(O) JdB x {0) JBi(O) 

It follows that u = in Bi(0) and that $3.0(0+) = 0. 
Last, in the case $ x o(0+) = —00 we obtain that 



lira S(x°,r)/r 2 = lim f r~ 3 ~" [ 



u 2 dH"- 1 

°) 



Taking a subsequence 

_ u(x° + r m x) 
Wm(X) :— . . 

that converges weakly in L 2 (dBi(0)) to w and setting T m := S(x° ,r m )/r m 2 , we infer 
from the monotonicity formula Theorem 14. II for large m that 

/ |Vu r J 2 < $3;o(ro) + / 2max(it rm ,0) + 2 u 2 Tm dU n ~ x . 

JSi(o) -/bi(o) Vasi(o) 
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Division by yields 

/ |Vu; m | 2 < T-H x o(r ) + T m x f 2max( Wm ,0) + 2 [ w 2 m dK 1 ' 1 . 

JBi(O) ifli(O) JdBi(0) 

Since |Aw m | < 1 in Bi(0), it follows that 

f |Vw | 2 < 2 f wl dH n - x , (5.12) 

JBi(O) JdBi{0) 

that u'o(O) = |Vu>o(0)| = and that wo is harmonic in £?i(0). From Lemma R~2l we infer 
that wo is a homogeneous harmonic polynomial of degree 2. 

Lemma 5.2 /n two dimensions, the only solution of that is homogeneous of degree 
2, is the trivial solution 0. 

Proof: Let u be a solution of (|1.2|1 on R 2 that is homogeneous of degree 2. Passing to 
the ODE y" + 4y = — X{y>o}> eacn component of {u < 0} must be a cone of opening 
7r/2 and each component of {u > 0} must for some r £ (0, +oo) be a cone of opening 
|arcsin(l/(4r))| < ir/2, a contradiction. 

Remark 5.3 Let u be a solution in a neighborhood of x° £ fl. If the point x° is non- 
degenerate of second order, then all points in some open neighborhood of x° are non- 
degenerate of second order. 

Proof: (0+) < implies by upper semicontinuity of 

x i * <M0+) 

that every point x £ {u = 0} in some open neighborhood of x° satisfies < J > a: (0+) < and 
is therefore non-degenerate of second order. 

6 Partial Regularity 

A result related to the following Corollary has been independently obtained in ^] The- 
orem 1.1]. 

Corollary 6.1 (Partial regularity in two dimensions) Let n — 2 and let u be a 
solution of ifi.il)) in fl that is non-degenerate of second order. Then for each K £C fl, 
the singular set K PI {u = 0} n {Vm = 0} contains at most finitely many points. 

Proof: Suppose this is not true. Then there is a sequence £1 R {u = 0} n {Vm = 0} 3 
x ,n — > x Q £ n {u — 0} n {Vu = 0}. Take a blow-up limit uq with respect to the fixed 
center x° such that dBi(0) contains a point of {uo = 0} fl {Vuo = 0}. By Proposition 
15.11 and Lemma 15.21 we know that uq is a homogeneous harmonic polynomial of degree 
2. This is a contradiction, since for a homogeneous harmonic polynomial of degree 2 in 
two dimensions the set {uq = 0} n {Vuq = 0} = {0}. 
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Lemma 6.2 Let u be a solution of in Q that is non- degenerate of second order, let 
x° G fl (~i {u = 0} R {V« = 0}, and let uq be a blow-up limit of 

_ u(x° + r m x) 

U m [X) '■ CV \ 

»j (x , r m ) 

in sense of Proposition I5.il Then for each compact set K C R n and each open set 
U D K n So there exists mo < oo such that S m R K C U for m > mo ; here So := {uo = 
0} n {Vm = 0} and S m := {w m = 0} n {Vit ro = 0} . 

Proof: Suppose towards a contradiction that S m H (if — [7) 3 £ m ->iasm^oo. Then 
5 G {mo = 0} n {Vuo = 0} n (K — U) , contradicting the assumption U D K Pi So. 

Proposition 6.3 (Partial regularity in higher dimensions) Let u be a solution of 
fl.ty) in Q that is non-degenerate of second order. Then the Hausdorff dimension of the 
set S — Q l~l {u = 0} n {Vu = 0} is less than or equal to n — 2 . 

Proof: Suppose that s > n — 2 and that H S (S) > . Then we may use ^] Proposition 
11.3], Lemma ffi. 21 as well as ^] Lemma 11.5] at Tt s -a..e. point of S to obtain a blow-up 
limit uo with the properties mentioned in Proposition 15.11 satisfying H S:Ca (So) > for 
So := {uq = 0} fl {Vuq = 0}. According to Proposition 15 . II there are two possibilities: 

1) (0+) = — oo and uo is a homogeneous harmonic polynomial of degree 2. But for 
such a polynomial H n ~ 2 ({u = 0} n {Vu = 0}) < +oo and we obtain a contradiction. 
Thus the second possibility has to apply: 

2) for some a G (0, +oo), auo is a solution of (|1.2fl on R™ that is homogeneous of degree 
2. In this case we proceed with the dimension reduction: By |14l Lemma 11.2] we find 
a point x G So — {0} at which the density in |141 Proposition 11.3] is estimated from 
below. Now each blow-up limit uqq with respect to x (and with respect to a subsequence 
m — ► oo such that the limit superior in 14, Proposition 11.3] becomes a limit) again 
satisfies the properties of Proposition 15.11 in addition, we obtain from the homogeneity 
of Mo as in Lemma 3.1 of |23j that the rotated Moo is constant in the direction of the n-th 
unit vector. Defining u as the restriction of this rotated Moo to R" _1 , it follows therefore 
that Ti s ~ l {{u = 0} n {Vm = 0}) > . 

Repeating the whole procedure n — 2 times we obtain a nontrivial homogeneous solution 
u* of degree 2 in R, satisfying 7i s ~(™ -2 ) ({m* = 0} n {Vm* = 0}) > , a contradiction. 

7 Lipschitz arcs in two dimensions 

In this section we show that the zero-set of the minimal solution consists in every second- 
order non-degenerate part of Q of finitely many Lipschitz arcs which end - if so at all - 
in quadruple junctions, meeting at right angles. 
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In order to do the analysis, we have to prove uniform Lipschitz regularity close to singular 
points. The difficulty is that convergence to the blow-up limit is not uniform at singular 
points. In I19j we used a novel intersection-comparison approach to obtain that close to 
the singular point the free boundary is uniformly the union of two graphs. In our case 
it turns out that the classical intersection-comparison method (also called zero-number 
technique or lap- number technique) is sufficient, when combined with a very elementary 
implicit function theorem argument. The proof of the following theorem is inspired by 

Theorem 7.1 (Unique blow-up limit) Let n — 2, let u be the minimal solution of 
in Q and suppose that u is non- degenerate of second order at x° G Q n {u = 

0} n {Vu = 0}. Then, as < r -> 0, and S(x°, r) := (r 1 -™ f gBr{x0] u 2 dH n -^ 2 , 

u(x° + rx) 
u r [x) := „ t 

o(x u , r) 

converges to p where (po U)(x) = (x\ — x"\)/\\x\ — x\ [|z> 2 (9Bi(o)) f or some rotation U. 
Proof: First, by Proposition [O] for any e > there is p > such that 

dist(u r , M g ) < e for r < p ; 

here M* := {{x\ - x 2 2 )/\\x\ ~ x 2 2 \\ L * {dBim }, M g := {q : q o V(x) = {x\ - x 2 2 )/\\x\ - 
^III^CaBiCo)) for some rotation V}, and dist(w r ,M s ) := mf qe M g sup xeBl ( ) \u r (x)-q(x)\. 
Denote by Ug the counterclockwise rotation of positive angle 8. If the statement of the 
theorem does not hold, then - by uniform continuity of t i — * n exp (_ t ) - there exists a 
sequence r rn { and rotations Ug 1 and Ug 2 satisfying \6\ — 6 2 \ = e S (0, 7r) as well as 

dist(w r2m o Ug 1 ,M*) < e and dist(u r2m+1 o Ug 2 ,M*) < e 

for m = 0, 1, 2, . . . 

Note that we may assume 9i > Q 2 . Now let U = U ej+s^ , oj = (#i — ^)/2 £ (0, 7r/2) and 

2 

define 

,, „, w(a;° + rU (cos 0, sine)) 

S(x^r) ' 

For each < r < r , the function <f>(r, ■) defines a function on the unit circle [— 7r, it). 
Inspired by applications of the Aleksandrov reflection (see for example J2]> P^| - 0) we 
consider now 

t(r,-6) :=#r,0)-#r,-0). 

Observe that £(r, 0) = ^(r, tt) = 0. In what follows we will prove that ||(r2 m ,0) < 
and ^|(r2, n +i,0) > for large m. The comparison principle (applied to the minimal 
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solutions S(x° ,r)<f>(r, 6) and S(x° ,r)<fi(r, —8) in the two-dimensional domain [0, ro)x (0, 7r) 
with respect to the original coordinates x\ and x%), tells us now that the connected 
component of {£ < 0} touching (r 2m ,0) intersects {r } x (0, n), and that the connected 
component of {£ > 0} touching (r 2m +i,0) intersects {ro} X (0,7r). It follows that {r } x 
(0, 7r) contains infinitely many connected components of {£ > 0} and {£ < 0}. On the 
other hand we know that, provided that ro has been chosen small enough, u ro o Ug 1 is 
close to {x\ — x\)/\\x\ — £2lU 2 (9Bi(o))> so €( r o,8) is close to ((9) := ci(cos(cj + 9) 2 — 
sin(cj + 6») 2 -cos(cj-6l) 2 +sin(cj-0) 2 ) = 2ci(cos(6» + lj) 2 - cos(6» - u) 2 ) in C x ([0,7r]). As 
the zeroes of £ are all non-degenerate, it is not possible that £(ro, ■) has infinitely many 
zeroes. 

Therefore, in order to finish the proof, we have to show that §|(r 2m ,0) < and 
§§(r 2m+ i,0) > for large m. As £(r 2m ,-) is close to C in C 1 ^,^]) and £(r 2m+ i,-) is 
close to — C in C 1 ([0, 7r]) we need only calculate C'(0) = — 8ci cos(cj) sin(w) < C2(w) < 0. 
So for e and ro sufficiently small (depending on e) we obtain a contradiction. 

Corollary 7.2 (Lipschitz arcs) Let n = 2, let u be the minimal solution of in 51 
and suppose that u is non-degenerate of second order at x° G 51 (~l {u — 0} D {Vu = 0}. 
Then {u — 0} consists in an open neighborhood of x° of four Lipschitz arcs meeting at 
right angles. 

Proof: By Theorem 17. II we know that for small S > and r 6 (0, S), the u r of Theorem 
17. II satisfies 

\d 2 (u r o CO | > ci in 52(0) n {|x 2 | > 1/8} 

and B 2 (Q) n {u r o U = 0} c Bi/i 6 ({N = M» ■ 

Thus (£ 2 (0) - S 1/8 (0)) n {w r o t/ = 0} is the union of four C^-graphs #t(j = 1,2,3,4) in 
the X2-direction, satisfying 

II Or)' llcg,-) < -rl|Viir[|c°(B a (o)) < C 2 , 
where Jj = [±1/2, ±1). Rescaling yields the statement of the corollary. 

8 Regularity of Local Minimizers 

Theorem 8.1 (Regularity of local minimizers) Let u be a minimizer of the energy 



in K := {v £ W 1,2 (B ro (x )) : v — u on dB ro (x )}. Then the free boundary d{u > 0} is 




locally in 51 an analytic surface and u 



{«>o} 



■, u 



{u<0} 



are locally in 51 analytic functions. 
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Remark 8.2 By Lemma \3. 7| this implies the same regularity for the maximal solution. 

Remark 8.3 The theorem also implies that local minimizers are locally in Q of class 
C ' . Usually regularity of the solution is proved before proving regularity of the free 
boundary, but here we do it the other way around. 

We start with some preliminary results. 



Lemma 8.4 Let n > 2, let u be a minimizer of the energy 



E(v) 



/ |Vu| 2 - 2max(v,0) 



in K :— {v £ W ' (B ro [x )) : v — u on dB ro (x )} , and suppose that Vu(i) ^ on 
{u = 0} n (B ro (x°)\ {x }). Then for w £ Cg°(B ro (x°)\B 5 {x )) we have 

< S 2 E{u) ■ H(w) = / |Vw| 2 - / TFTa™ 

Proof: We define 



B ro (x°) J{u=0}nB rQ (x<>) lVu| 



E e (u)= [ \\Vu\ 2 - le {u) 

JB ra (xO) 2 

where 7 e (w) is an approximation of max(u, 0) such that "f'J(u) = 1/e if u £ (0, e) and zero 
otherwise. We see that for w £ (B ro (x°)\Bs(x )) we have 

1 

7 

where 

If t 2 

A \ = 72 / t' Vw I 2 ~ (T«( u + tw ) - 7«(«) ~ HW)w) 

1 JB T0 {x°) 1 

can be rewritten as 



2 (E e (u + tw) - E e (u) - tSE e (u){w)) = A\ 



A\ = ( f^|Vw| 2 - [ da [ da rfj (u + atw) 

JB rn (x°) V 2 Jo JO 



w 2 



r (x°) 

l\Vw\ 2 - [ da [ - [ w 2 . 

B rQ {x a ) 2 Jq Jq £ J B ro (x°)n{0<u+atw<£} 



By the co-area formula, we obtain for small t that 



l 



A\-> A l = / \\Vw\ 2 - I da I da 



lB ro (x°)'2 Jo Jo J{u+atw=0}nB ro (x°) \^ ( u + Crtw)\ 

as e — > 0. We conclude that 

< 4 ( E ( u + tw ) - E ( u ) - tSE(u){w)) = Ai . 
t A 
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Last, we take the limit t — > and obtain 

1 JB ro (x°) 1 1 J {u=a}nB ro (x°) I Vw 



Lemma 8.5 If u is a solution in B ro (0) of satisfying Vu(0) = 0, then there exists a 
constant C < oo such that 

\Vu(x)\ < C|x|logl/|x| inB ro/2 (0). (8.13) 

Proof: From Chemin we infer that Au 6 L°° implies u £ C 2 . Thus Vu G C*, and 
there exists a constant C < oo such that 

|Vu(as)-Vu(y)| < C|x - y|(l + log l/|as - y|) in S ro/2 (0) 

which proves the Lemma. 
Proof of Theorem EH 
Step 1 (Dimension n = 2): 

From Lemma 13.61 and Corollary 16.11 we know that locally the free boundary is either 
a C 1 -arc - in which case the gradient is non-zero on the free boundary - or a cross 
composed of 4 Lipschitz arcs meeting at right angles - in which case the gradient is zero 
at the center x 1 of the cross. We want to show that such a cross is impossible for local 
minimizers. We may assume x 1 = 0. 

From Lemma \8. 41 and (|8.13() . we deduce that for some constant c\ > 

< f |Vu,| 2 -ex I | M * V for w e C o °°(B ri (0) - B 5n (0)) . 

JB ri (o) J{u=o}nB ri (o) F|logl/|a;| 

(8.14) 

We now consider Ws{x) — <p(x) — (j)(j) where cj) € C^(B ri (0)) such that = 1 on_Bn(0). 
It follows that J B ^ \ Vw$\ 2 < C2 < oo, and using for large i S N the regularity of 
{u = 0} (~1 (B 2 -i(0) — B 2 -i-i(0)) as well as the closeness to the rotated cross we obtain 

C 3 / 1 1 1 1 1 71 , wj > - I dr > +00 

y{«=0}ns n (0) |x|logl/|x| Js rlogr 

as 5 — > 0, a contradiction to the boundedness of 

Therefore the cross is not a local minimizer, and in dimension n — 2 the free boundary 
is locally in £1 a C^-arc for each local minimizer. 
Step 2 (Dimension n > 2): 
We proceed by induction. 

We assume that we have proved that for local minimizers, the free boundary is smooth 
up to the dimension n — 1 > 2. 
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Now we cannot have an accumulation of singularities in dimension n: Blowing up at a 
limit point and blowing up a second time at a singularity x ^ of the blow-up limit, we 
would obtain as in the proof of Proposition 16.31 (see also |231 Lemma 3.1, Lemma 3.2]) 
a local minimizer with a singularity in dimension n — 1. Thus singularities are isolated, 
and every blow-up limit of u at each singularity is a harmonic polynomial of degree 2 
whose gradient vanishes only at one point of the 0-level set. Thus we still have (|8.14() . 
and the free boundary is for large ieNon B 2 -i(0) — B 2 -i-i(0) close to the zero level 
set of a homogeneous harmonic polynomial Pi of degree 2 satisfying by Proposition 15.11 
\Pi(x)\ < C 4 \x\ 2 and |VPi| > c 5 |x| on (R™ - {0}) n {P t = 0} where C 4 < oo and c 5 > 
do not depend on i. 

Now we choose ws(x) = \x\~(^~\(f>(x) - 0(f)) where <j> £ Cg°(B ri (0)) and <j) = 1 on 
Bn(0). We obtain J B , Q » |Vu>«5| 2 < Ce, and using the regularity of {u = 0} fl {B 2 -i (0) — 
B 2 -i-i(0)) as well as the closeness to {Pi = 0} we see that 

C 7 f \ w 2 5 >- / 2 dr - J— finf f i dU n - 2 \ -> +oo 

as 6 — > 0, a contradiction to the boundedness of wg- 

Therefore the local minimizers have no singularities in dimension n, and the free boundary 

is locally in a C 1 -surface. 

Step 3 : Analyticity of the free boundary 

We obtain analyticity of the free boundary as well as analyticity of u \ ^ u>0 y u \ { u< o\ as 
in the proof of Theorem 4.1 in Chapter 6 of JHI- See also |151 Theorem 3.1']. 



9 The Cross Singularity 

The reader may wonder whether there exists an example of a singularity for the maximal 
solution (and thereby a counter-example to the W /2 '°°-regularity in this unstable prob- 
lem). We have at this moment no conclusive answer, but the following formal asymptotic 
expansion suggests that the cross may be a possibility: 

Lemma 9.1 (Formal asymptotics) 

Let us assume the existence of a solution u in -Bi(O) such that u satisfies for 

x\ = r cos a, x 2 — r sin a 

u(~xi,x 2 ) = u(xi,x 2 ) = u(xi, -x 2 ), 
and the free boundary in the set -E>i(0) n {x\ > 0, x 2 > 0} is given by 

a =7(1 + MO), P=—!—- 
4 — log r 
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Moreover we assume that 



u>0 in < a < -(1 + <j>{p)) 



u<0 in -(1 + (}>{p)) < a < tt/2 



Then formally <ft and u satisfy 



m = -^ + o(p 2 ) 



and u{x) = — (x\ - x\){- log \x\) + 0{\x\ 2 ) . 



Formal proof: 
We set 



u(xi,x 2 ) = < 



' r 2 (-log(r)) z+(p,0+) for 9+ = - 



-r 2 (— log(r)) z (p,6 ) for 



if -(i + ^x^^. 



By continuity of u and Vu on the free boundary, we have 

z+ = z~ = and z+ = ( 1 + ^f] ) zl for 9+ = 0" = n/4 . 

By the symmetries of u we also have 

z+ = z g = for 0+ = 0- = . 
Moreover z = z + satisfies for < 9 = 9 + < 7t/4: 

zee + 4(1 + Hp)) 2 z + (p- Apz) (1 + 4>{p)) 2 + p 2 I[z, 4>] = 

where 

I[z, 4>] = 4 ((1 + c/>(p)) 2 z p - ez e 4>'{p){l + 0(p))) 

+ p 2 ((1 + 0(p))% P + Ap)(^ + 20z e ) - (1 + cj>{p)) {29z p g4>'(p) + 9z e 4>"(p)}) 

and for < 9 = 9~ < 7r/4, the function z = z~ satisfies a similar equation with <\> 
replaced by —<j> and without the term p(l + 4>{p)) 2 , i.e. 

zee + 4(1 - 4>{p)) 2 z - Apz(l - 0(p)) 2 + p 2 I[z, -0] = . 

Let us now introduce the formal asymptotic expansion 

' 4>{p) = 4>° + P& 1 + p 2 <t> 2 + ... 

< 
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where 

<f)° = 0, z ±,o (0) = A°cos(26) and A > by assumption. 
For the order terms we obtain (for < 8 < 7r/4) 

4'° + 42 ± '° = ) ^°(0) = 0, z^°(V4) = 0, 4'°(7r/4) = V°(tt/4) , 

which is compatible with the assumptions. For the order 1 terms we obtain (for < 9 < 
vr/4) 

zjg 1 + AZ+' 1 + (80 1 - 4)z+>° + 1 = 0, Zgg 1 + 4z ~ > 1 + (-80 1 - 4)z+>° = 0, 

[ zp'iO) = 0, z±> l (w/4) = 0, z+'V/^ = ^' 1 (tt/4) - 2^z e >°(n/4) . 
Thus 

z+^tO) = c,os(20) 4- (1 - 26 1 )A° 9sm(20) - 



z +A (6) = A+' 1 cos(26>) + (1 - 2</> 1 )A 6>sin(26>) - \ 
and z-'\0) = A-- 1 cos(26) + (1 + 2<t > 1 )A°6 sin(26>) 



with a!) 1 = -i A = ^- and - A - ' 1 = A. 

~ Z 3 Z7T 7T 



10 Open Questions 

The most urgent remaining questions in this context are, whether the cross singularity 
can be proven to exist, and whether there are examples of second order degeneracy. 
Another interesting point is whether methods similar to those used in this paper (possibly 
combined with arguments as in pQ) can be used to prove regularity in the composite 
membrane problem (see [H], ^UJi OH)- 
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